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Abstract: In spacetimes with compact dimensions there exist several black object solu- 
tions including the black-hole and the black-string. These solutions may become unstable 
depending on their relative size and the relevant length scale set by the compact dimen- 
sions. The transition between these solutions raises puzzles and addresses fundamental 
questions such as topology change, uniquenesses and cosmic censorship. Here, we consider 
black strings wrapped over the compact circle of a d-dimensional cylindrical spacetime. We 
construct static perturbative non-uniform string solutions around the instability point of 
a uniform string. First we compute the instability mass for a large range of dimensions, d, 
and find that it follows essentially an exponential law 7^^, where 7 is a constant. Then we 
determine that there is a critical dimension, d^: = 13, such that for d < d^ the phase tran- 
sition between the uniform and the non-uniform strings is of first order, while for d > d^, 
it is, surprisingly, of higher order. 



In 4d the static uncharged black hole(BH) solutions with a given mass are stable 
and unique. However the fundamental theory of nature, which as now believed by many, 
is the string/M-theory contains more than four dimensions. In this situation the phase 
space of massive solutions of General Relativity is much more rich and varied. Several 
phases of solutions exist and transitions between them may occur. For concreteness, we 
consider the background with a single compact dimension, i.e. with the topology of a 
cylinder, M*^"^'^ x S^. The coordinate along the compact direction is denoted by z and its 
asymptotic length is L. The problem is characterized by a single dimensionless parameter 

/i := GrfM/L'^-^, (1) 

where Gd is the d-dimensional gravitational constant and M is the mass. 

Gregory and Laflamme (GL) |2| discovered that the uniform black string (i.e. a 
d — \ Schwarzschild solution times a circle, which is the large mass solution) develops a 
dynamical instability if the compactification radius is "too large". Their interpretation 
was that the string decays to a single localized BH. In this case the horizon pinches off 
and the central singularity becomes "naked". By now there is a rapidly growing amount 
of the literature on the subject [3-21]. In particular the scenario of GL was questioned by 
Horowitz and Maeda (HM) who, on grounds of the classical "no tear" property of the 
horizons, argued that horizon pinching is impossible and hence a decaying string settles to 
another stable phase - a non-uniform black string(NUBS). However, a (partial) evidence 
against that has come from Gubser Q who in 5d studied perturbative NUBSs emerging 
from the GL point. He showed that such solutions are too massive and have too low an 
entropy to serve as an end-state of a decaying critical string. Namely, the transition to 
this NUBS is of first order and it is again unclear what state is accessed by the classically 
decaying GL string. Wiseman reached the same conclusion by constructing the NUBS 
solutions numerically in 6d in a fully non-linear regime^. However, in this paper we discover 
that the transition to NUBS can be smooth depending on spacetime dimension. 

Generalizing Gubser 's 5d procedure |Q, which is a version of the "marginal stabil- 
ity" method, we construct numerically d-dimensional static perturbative NUBS solutions 
around the GL point. First, we note that the GL instability mass exhibits to a good ac- 
curacy an exponential scaling with d. Moreover, we find that there is a critical dimension, 
dif = 13, below which the uniform-nonuniform strings transition is of first order. I.e. it 
is qualitatively similar to what Gubser has found in 5d. However, above d^: the NUBS 
solutions emerging from the instability point have a lower mass and a larger entropy than 
those of the critical string. Namely, the transition between the phases can be continuous^. 

Hence the NUBS state is accessible by an unstable uniform string. In this case, the 
horizon would not pinch off at the GL point. Our result suggests, however, that the horizon 
fragmentation during the classical decay can be avoided only for d > d^. This is a rather 
curious development since the original HM argument was dimension independent. It should 

^In 5d, jl^, |l^ could be regarded as additional circumstantial evidences contra the HM claim. 
^This is consistent with the prediction of a critical dimension d = 10 at the "merger point" of this 
system, where the string and the BH branches merge M. 
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be noted, however, that the central issue of whether any unstable string must decay to a 
string remains unresolved even for d > d^,. 

The most general ansatz for static black string solutions is 

/ = 1 - l/r''-\ (2) 

where A, B and C depend on r, z only. When these functions vanish the metric becomes 
that of a static uniform black string with the horizon located at r/^ = 1. 

Gubser Q has considered static NUBS solutions that differ only perturbatively from a 
uniform black string. Since the method was described in detail in the original paper |Q] and 
then in Q we mention only the most important points. Gubser developed a perturbation 
theory considering the expansion of the metric functions in powers of A. This A parametrizes 
the NUBS branch that joins the GL point in the limit A ^ 0. The expansion has the form'^ 

oo 

X = ^\''Xn{r)cos{nKz), 

n=0 

oo oo 

Xn{r) = J2 A'^^„,p(r), K = Y^ (3) 

p=0 q=0 

for X = A,B,C with Xo,o = 0; and K = 21: jL. 

Upon substituting (^) into the Einstein equations, R^^_l, = 0, a finite set of ODEs is 
generated at each order of the expansion^ . Gubser's method is very accurate up to the 
third order in A. Following the original procedure we restrict our computations up to 
O(A^). Nevertheless, interesting results are already obtained here. Actually, the third 
order is precisely what one needs to determine the smoothness of a phase transition. 

As discussed in Q the perturbation theory contains a "scheme" dependence that seems 
to correspond to different parameterizations of the non-uniform branch. Originally, fixing of 
the "scheme" was achieved by fixing the constants c„^p := Cn,p{rh)- Still, other "schemes" 
can be used. For example in [Q] the asymptotic length of the compact circle was held fixed, 
K = const, but the constants c„^p were allowed to vary. In fact different "schemes" all 
produce the same scheme-independent results, like e.g. the dimensionless mass (111). Here 
we choose to work in the "standard scheme", as it is referred in Q, by fixing Cn,p = for 
n > 1 and ci^ = 1. 

Once the metric functions are known, various thermodynamical variables can be com- 
puted. Asymptotically the spacetime (|2|) is characterized by two charges [16, 19 1 - the mass 



and the tension of the black string. By making a Kaluza-Klein reduction in the z direction, 
Xn,p in (^) are observed to be massive modes for n > and they are massless otherwise. 



^Gubser used the "non-uniformity" parameter, A ~ l/2(i?max/i?min — 1) where -Rmax and -Rmin refer 
to the z-dependent Schwarzschild radius of the horizon. Hence A coincides with this A only at the leading 
order. It was shown subsequently in [ p^ , pj| that a good order parameter that allows to put black strings 
and holes on the same phase diagram is not A, which is undefined for the latter, but the scalar charge of 
the dilatonic field. However, for our current purposes A may be left unspecified. 

"^See e.g. for derivation of the Einstein equations in a very similar case. 
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Only the latter contribute to the asymptotic charges since the former decay exponentially. 
Up to O(A^) the relevant massless modes are Xo,i- Asymptotically, they fall off as inverse 
powers of r. We denote the coefficients of the leading terms by X^o- It is convenient to 
define the variation of the charges of a non-uniform string with respect to a uniform one. 
According to at the leading order these variations read^ 

5M/M = -2 [A^ + B^/{d - 3) ] A^, 

5T/T = -2 [A^ + {d- 3)Soo] A'. (4) 

We also compute the variation in the temperature, 5T /T = exp[^ — B] — 1, and in the 
entropy, 5S/ S = exp[i? — [d — 3)C] — 1, which are evaluated at r = 1. 

Finally, defining the variation of K, 5K/K := {ki/ko))?, we determine the dimension- 
less, scheme-independent variables by multiplying the dimensional quantities by suitable 
powers of K. By doing so we obtain for our variables 

5ii/^i = 6M/M + {d- A)5K/K := 7]iX^ + . . . , 
5t/t = 5T/T + {d- A)5K/K := nA^ + . . . , 
89/9 = 6T/T - 5K/K := ^lA^ + . . . , 

5s /s = 6S/S + (d - 3)5K/K := siA^ + . . . . (5) 

Incorporating the first law as in we evaluate the entropy difference between the 
non-uniform and uniform strings with the same mass 

'S'non— uniform . \ 2 , \ 4 , 

^; = i + CJiA + <T2A +..., 

'-'uniform 

d-4 d-3 1 \ 

= - d^3'' ' = -w^) r ^ d^^""' J ^ ^ 

The vanishing of ai is ensured by the first law at the leading order (where L = const) Q. 
We verified that to a good ( ^ 1%) accuracy, ai ^ for our solutions. Thus, the entropy 
difference (|6|) arises only at O(A^). 

At each order of A we solved the ODEs numerically[^2|. We were able to exactly 
reproduce the numbers found so far in the literature: for 5d in |^ and for 6d in An 
indication of the accuracy of the method is gained by varying the "scheme" [Q] , by altering 
co,i = 0, ±1. The resulting variation in (^J^ gives an idea of the numerical uncertainty. 
For small d's the accuracy of our calculation is high, being about 0.5% in r/i and 1% in a2- 
For larger d's the method is somewhat less accurate, yielding 5% and 6% variations in rji 
and (72 respectively, for d = 16. This has to do with the steep asymptotic fall off of A and 
B in which the leading terms decay as r~^'^~^\ while C falls off only as 1/r (in 5d the fall 
off is log(r)/r). Hence, the accuracy in extracting the coefficients Aoo,Boo, that contribute 
to r]i and o"2, decreases for large d. 

The critical mass. The calculation in the linear order in A yields the mass of the critical 
string, since the leading order of (0) corresponds to the static GL mode. We performed the 



^We use units in which Gjv ~ Gd/L = 1. 
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Figure 1: The relative difference between tlie mass and the fit (^), 0.477'^, as a function of d. 
For lie this difference is zero with the spread of about 0.8% magnitude, giving approximately 2.1% 
variations in /ic itself. 



calculations in d = 5, . . . , 16, 20, 30 and 50. For d < 10 we confirm a very good agreement 
with the original GL resultsQ, presented in their FIG. 1. Note, however, that the methods 
are very different. For the entire range of d we find that the critical mass is remarkably 
well approximated by 

II, oc Y (7) 

with 7 ~ 0.686 being a constant, and the prefactor is approximately 0.47, for the specific 
definition of mass (|l|). In FIG. [l| we plot the relative difference between the logarithm of 
the critical mass and the fit (0). It is clearly seen that log(/ic) is linear for all d. There is 
still room for a weak d-dependence, of order 2.1%, around the dominant scaling (Q). We, 
however, could not extract this residual dependence. 

To get an insight into this behavior (|7|) we compute the mass of a uniform black 
string whose entropy is equal to that of a single BH with the same mass. First we 
compare the entropy of the black string with that of a d-dimensional Schwarzschild BH. 
The corresponding entropies read = Aj__2/{4:Gd) and S'sstr = ^d-3-^^/(4Grf) where 
Ad := J7d[167rM/(dOrf)]'^/('^-i) and fla is the surface area of a unit S"^ sphere. Equating 
these, S^^{iJ,) = SBStiil^'), we solve for the mass 

„(o) ^ 1 ^r^ (rf-3)('^-3)(-^-3) 

^ 167rl^j:4(^_2)(<i-2)(d-4)- W 
Actually, we can do slightly better by using the analytical formula for the entropy of 
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Figure 2: The trends in tlie mass, A*non-uniform/Muniform := 1 + f]iA^ + . . ., and the entropy, 
•Snon- uniform /'S'uniform := 1 + (^2^'^ + • • ., shifts between Uniform and non-imiform black strings. The 
key result is the sign change of rji and CT2 above = 13. 



small BHs on cylinders derived recently in |21 



_ c(0) 

'-'bh ~ '-'bh 



- C(^-3)16^^ 2 



(9) 



where C,{x) is Riemann's zeta-function. This formula reflects the leading order corrections 
to the Schwarzschild metric due to compactification^ . It implies that for a given mass 
the entropy of a "caged black hole" (a BH in a compactified spacetime) is larger than the 
entropy of a Schwarzschild BH. The mass /x^^^ corresponding to equality of the entropies 
is then obtained by solving the equation S^^{ii) = S■Q■s,\.v{^J) ■ 

We add to FIG. || the plots for these masses. In contrast to log(/ic) the logarithms 
of ^^^^ and /i^^^ have a non-linear dependence on the dimension for small d's. They do, 
however, become linear (with a different slope) for d S> 10. Here we already see a hint of 
a critical dimension - looking at the difference between /ic and its estimator (either ix^'^^ or 
/i^^)) one notices a change of sign at about d ^ 12.5. This suggest that for d ^ 13 the BH 
state is entropically favorable over the string state only for ^ < /ic- 

From a sudden to a smooth phase transition. Performing the computation in higher 
orders, up to 0{\^), we obtain the variation in the variables (^) and entropy The 
results for 771 and (T2 are depicted in FIG. ^. One observes that rji is initially positive for 
d = 5, reaches a maximum at d = 10, and becomes negative for d > 13. Then it continues 
to decrease and in fact it drops increasingly faster with d, as indicated by the growing 
distances between subsequent points in the graph. The pattern for cj2 is similar but with 



^The perturbation theory is constructed in powers of /i <C 1. 
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the opposite sign^. 

The key phenomena is the appearance of a critical dimension, = 13, above which the 
perturbative non-uniform strings are less massive than the marginal GL string. Moreover, 
their entropy is larger than the entropy of the uniform string with the same mass. It is 
important that ryi and (T2 change signs simultaneously. 

As for the other variables, we find that the trend in the entropy shift, si, is qualitatively 
similar to the behavior of 7?i - it is positive for d < and it becomes negative above d*. 
For the variation of the temperature we note that below d^, the NUBS is "cooler" than 
the uniform one and above c?* it is "hotter" . We find that the tension of the non-uniform 
strings is lower than that of uniform ones. This is in tune with the expectation that the 
uniform black string has a maximal tension, and that the tension vanishes for small black 
holes [|^, 23 1 . In addition, we observe the ratios r/i/ri and r/i/si to be discontinuous near 
d*. Note also that in FIG. ^ we plot the coefficients of the mass and the entropy shifts. To 
obtain the physical variations these and others coefficients must be multiplied by suitable 
powers of A. 

To summarize. While we have found the dependence of the critical mass on the di- 
mension we do not have at present an explanation for the scaling We believe, it gives 
us some insight into the nature of the GL instability and it probably is connected with the 
thermodynamical instability of the system|14|. However, it is the appearance of a critical 
dimension, d^,, that can perhaps be regarded as our main result. It implies that above 
the critical string can smoothly evolve into the NUBS phase. For d < d^ the transition 
between the two phases is of first order. 

The continuous transition above d^, suggests that the NUBS phase can be a natural end 
state of the GL instability. Indeed, a uniform string loosing its mass by evaporation and 
encountering the instability at /Uc can smoothly evolve to the non-uniform state keeping its 
singularity covered by the horizon. Already from FIG. |^ it could be inferred that above 
d ^ 13 there can be a branch of solutions between the uniform strings and the BHs. We 
believe that the NUBS state is a reasonable candidate for this "missing link". 

As the mass is further radiated away two scenarios may be proposed: (1) The NUBS 
branch extends to an arbitrary small mass. A black string evolves along this branch 
probably increasing its non-uniformity all the way down to zero mass. In this case the 
cosmic censorship would be held(at least until the final stages of evaporation); (2) A NUBS 
becomes unstable at a finite mass where the horizon fragments and a localized BH forms. 
This may lead to a compromize of the cosmic censorship, much like in the d < d^ case but 
for a mass smaller than ^c- The transition between a NUBS and a BH can be sudden or 
smooth depending on the relative values of the instability masses for these states. Note 
that a NUBS branch that extends to zero mass or becomes unstable even earlier on a phase 
diagram is conceptually the same. The main difference is whether the naked singularity 
shows up before the end of evaporation or not. 

To address these intriguing issues it would be a very interesting future task to construct 
in a fully non-linear regime, like in Q, the branch of NUBSs that we found here. In 

''in fact we also did the computation in d = 20 finding the same trends. However, the numerical errors 
were of order 20% so we regard this case as an indicative only. 
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particular, it is interesting to determine for how low a mass this branch drops, would the 
horizon try to pinch off forming a cone-like "waist" |9| and whether the topology tends 



to change. In addition we expect that a time evolution of the critical string, like in |12|, 
should confirm a nice decay for d > d^^. 

In this work we have considered black strings in a cylindrical spacetime. We believe 
that the critical dimension phenomena is general and will hold for more general backgrounds 
with additional compact dimensions even if the specific value = 13 would change. 

I thank B. Kol and T. Piran for stimulating discussions and valuable remarks on the 
manuscript and Mu-In Park for pointing out an error in the /i^^^ calculation. 
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